BKT-like transition in the Potts model on an inhomogeneous annealed network 
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We solve the ferromagnetic g-state Potts model on an inhomogeneous annealed network which 
mimics a random recursive graph. We find that this system has the inverted Berezinskii-Kosterlitz- 
Thouless (BKT) phase transition for any q > 1, including the values q > 3, where the Potts model 
normally shows a first order phase transition. We obtain the temperature dependences of the order 
parameter, specific heat, and susceptibility demonstrating features typical for the BKT transition. 
We show that in the entire normal phase, both the distribution of a linear response to an applied 
local field and the distribution of spin-spin correlations have a critical, i.e. power-law, form. 
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I. INTRODUCTION 

Among numerous unusual phenomena which were ob- 
served in complex networks 0, H, H, 0, [f| , one effect is 
astonishing. In a wide range of growing networks, a phase 
transition of the birth of the giant connected component 
demonstrates the Berezinskii-Kosterlitz-Thouless (BKT) 
singularity @, 0, i, i, M, d d d d , i.e., the size S 
of the giant connected component has the critical singu- 
larity S oc exp (—const /V8). Here the appropriate pa- 
rameter S characterizes the deviation from the critical 
point. Recall that in cooperative models on lattices, an 
(infinite order) BKT transition is possible in very specific 
situations (e.g., in the two-dimensional XY model) where 
critical fluctuations are quite strong [ll| [l6| , on the lower 
critical dimension. In contrast, networks are extremely 
compact, infinite-dimensional objects (small worlds) and 
so can be described by, in essence, mean field theories. 
Consequently the reason for the BKT singularity in grow- 
ing network systems strongly differs from that for the 
BKT transition in solid state systems. It was argued that 
it is the inhomogeneous architecture of some growing net- 
works that results in the BKT singularities. Note that 
in equilibrium networks with any degree distributions, 
where all vertices are statisticall y equivalent, this transi- 
tion is impossible [13, El El Hfl Hit H [H, M, EE El • 

Interestingly, a similar critical behaviour was earlier 
(already in 1990!) found, proved, and analysed in the 
ferromagnetic Ising model for a chain of spins with a 
specific large-scale inhomogeneity of long-range interac- 
tions [13, [2i|. (In Ref. 28], the BKT singularity was 
also revealed in another cooperative model on the same 
substrate, where, instead of the Ising spins, there were 
continuous scalar fields.) In a network context, actu- 
ally the same system was studied in Ref. [2?|, where a 
random recursive graph (i.e., with quenched disorder) 
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was approximately substituted by a graph with annealed 
disorder — an annealed network. (By definition, in the 
random recursive graph, at each time step, a new ver- 
tex is attached to a randomly chosen existing vertex.) 
The resulting model may be analysed exactly, and it is 
equivalent to the Ising model on a deterministic fully con- 
nected graph with inhomogeneous Ising interactions be- 
tween spins (see Fig. [I}. The BKT singularity in the Ising 
model on another, hierarchically organized, deterministic 
network were found in Ref. [30] . 

One should stress that on inhomogeneous networks, 
the BKT transition was found in the cooperative models 
(the Ising model, percolation) which have a second order 
phase transition in high-dimensional lattices and the clas- 
sical random graphs. In contrast, in the present paper 
we solve the g-state ferromagnetic Potts model, which 
has a first order phase transition on the classical random 
graphs if q > 3 [311 ]. That is, there is a jump of the order 
parameter and hysteresis if the substrate of the model 
is an infinite-dimensional lattice. We show that even in 
this, quite different situation, a BKT-like infinite-order 
phase transition is realised on the annealed variation of 
the random recursive graph, and the jump and hystere- 
sis are absent. We calculate a set of observables and 
find typical BKT singularities. Furthermore, we obtain 
distribution of the linear response to local field — the dis- 
tribution of "correlation volume" — and find that it has a 
power law form in the entire normal phase. Recall that 
in cooperative models on equilibrium networks, e.g., on 
uncorrelated networks, this correlation characteristic is a 
power law only at a critical point. We also obtain the 
distribution of pair spin-spin correlations, which was re- 
cently introduced in Ref. [321 ] , and again find a power law 
in the entire normal phase. Thus we generalize results 
obtained for the Ising model on this network to the case 
of the Potts model with an arbitrary number of states. 
The revealed critical behavior dramatically differs from 
that of the Potts model on the configuration model [24j |. 
Technically, our derivation directly follows Ref. [29] ]. 

This paper is organized as follows. In Sec. [Til we intro- 
duce the model, in Sec. IIII[ we present our main results, 
and in the following sections derive and discuss them in 
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FIG. 1: The deterministic fully connected graph which is 
equivalent to the asymmetric annealed network under study. 
Here t = 3. The values of the Potts couplings are shown on 
the edges. 

detail. 



II. THE MODEL 

The Hamiltonian of the ferromagnetic g-state Potts 
model is 

t 

« = 5^(^,1), (1) 

(ij) i=0 

where the coupling is set equal to 1, (ij) indicates sum- 
ming over the coupled spins, and Si = 1, 2, q. S(si, Sj) 
is the Kronecker delta symbol, and Hi is the external 
field. At q = 1, this model is equivalent to the site per- 
colation problem [33|, and at q = 2, it is the Ising model. 
For these two values of q, at the standard mean-field 
regime, the Potts model have a second order phase tran- 
sition. For q > 3, the standard mean-field theory gives 
the first order phase transition. 

We place the Potts model on the following asymmetric 
annealed network. Vertices are labelled i = 0, 1, 2, . . . , t, 
as in a growing network. Each vertex, say vertex i, 
have a single link of unit strength to previous vertices 
j = 0, 1, 2, . . . , i — 1. The second end of this connection 
frequently hops at random among j = 0,1,2, ... ,i — 1, 
which just means specific asymmetric annealing. The 
resulting network is equivalent to the fully connected 
graph with a large-scale inhomogeneity of the coupling 
(see Fig. [T]). The mean-field theory for this network is 
exact, and it also mimics the more difficult case of the 
random recursive graph — the case of a quenched disorder 

In other words, in this model, the spin which was born 
at time i, has equal coupling 1/i to each of the older 
spins. That is, the Hamiltonian of our system is 

n = - J2 %^-X><^,i)- (2) 
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FIG. 2: The order parameter M (a), the specific heat C 
(b), and the susceptibility in zero field, x> (°) versus inverse 
temperature (3 = 1/T at different q in the g-state Potts model 
on the asymmetric annealed network. The numbers on the 
curves indicate corresponding q. The critical point is f3 c = 
1/Tc = <?/4. x{Pc-) = 1, x(/3c+) = 3. 



III. MAIN RESULTS 



Here we briefly list our results obtained for the infinite 
annealed network t — > oo. We demonstrate that for the 
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all g > 1, our system has an infinite order phase transi- 
tion at the critical temperature 



(3) 



IV. MEAN-FIELD EQUATIONS 



It is convenient to rewrite the Hamiltonian @ as 



In the following we use = 1/T for the sake of con- 
venience. The resulting mean spontaneous "magnetiza- 
tion" [M = q(S(s,l) — l/q)T/(q— 1), where (-)t means 
the thermodynamic averaging] is of the typical BKT form 



M{0) 



q 



(4) 



near the critical temperature [see Fig. [Ha)]. Here 
c = 1.554.... Note the difference by factor 2 of the 
value T c (q=2) from the corresponding result for the Ising 
model. This usual difference (see, e.g., Refs. 0, HH) 
disappears after the proper rescaling of the coupling con- 
stant in the Hamiltonian ((T|) to arrive at the Ising Hamil- 
tonian with the coupling constant 1/j. 

The specific heat Ch=o{T) = C(T) per spin is zero in 
the normal phase, and, below the critical temperature is 



C (0) = ~J28~~ -77T3 ex l» 



(/?-g/4) 3 



-2tt/ 



V 



(5) 



where (7rce) 2 /128 = 1.37. . . [see Fig.^b)]. 

We also find the susceptibility of the system in zero 
field: 

X (0 > g/4) = q/0-l if - g/4 « g/4, 

X (/3 < g/4) = (l - y/1 - 40/q) / (l + a/1 - 4/3/?) • (6) 

There is a finite jump of the susceptibility at the phase 
transition point: x[/3=(g/4)-]=l and x[/?=(g/4)+]=3 
[see Fig. He)]. 

The distribution of the linear response Xi to a local 
applied field (i.e., the correlation volume distribution), 

Xi) ) /t, is a power law in the entire 



normal phase: 



P(^)ocm" [1+2/(1 -^ 4/3/<?)1 - 



(7) 



Below the critical temperature, this characteristic of cor- 
relations rapidly decreases. Thus there is a contact of the 
"critical" phase with the power law decreasing distribu- 
tion and the phase with a fast decrease of P (fi). At the 
phase transition point, P (/x, — g/4) oc /i -3 . 

Furthermore, we calculate another distribution func- 
tion P{v) = (2/t 2 ) (j2l,j 5 ( v ~ Xij)), that is the distri- 
bution of correlations ~ Xij m the neighboring spin pairs 
O (Xi — X) j Xij) > an d above the phase transition, we 
find a power law with the same exponent: 



(8) 



1 q * 



0<i<J<t P=l 



i=0 



Let us assume small fluctuations of 5(si,p) from their 
mean- field values h ip . In the next section we shall show 
that this usual (mean-field theory) assumption leads to 
exact results. So, one can substitute 



5(si,p)5(sj,p) — ► -hiphjp + hi P 5(sj,p) + h JP S(si,p) 

(10) 

into the Hamiltonian @, where hi p — (6(si,p)}. This 
results in a linear effective mean-field Hamiltonian. The 
partition function of this Hamiltonian is: 



Z = exp E 



t q 



IlE ex p 



p=l 0<i<j<t 
' ' i-1 



/?(E — + .t ^) + 0HiS(p,l) 



=0 1 i=»+i 3 



(11) 



For large t, one may pass to the continuum limit: hi P 
h p (x — i/t). Then the partition function is 



l l 



Z = exp I —0 / / dxdy 



P =i 



h p (x) h p (y) 



x 



t q 

nE ex p 

!=0p=l 



3\ljdyh p (y) + Jdy f ^-}+ .UKrUnp.l) 

x 



(12) 



and therefore the free energy F = —0 1 In Z is 



l l 



F = tJ2 dxdy 



hp (x) hp (y) 



x 



jj Jdx In } ^2 ex P ( ~ J d V h p (y) + J 



dy hAy) 



0H(x)6(p,l) 



(13) 
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which leads to the self-consistent equations for the mean 
values h ip = J2{ s =1 2...0} 5(s l ,p)e l3n{s ^ / Z t : 



hi (x) 
1 



A(x) 
h p (x 
1 



exp 



i[- [ X dyhi(y)+[ dy^- ) -M7/(.ri 

x J o Jx y 



A(x) 
where 
A(x) 

exp 



exp 



P[- dyh p (y)+ dy 



,p = 2,...q, (14) 



P 



dyhi{y) + / dy 



hi(y) 



f3H{x) 



^exp 

p=2 



P[~ dyh p (y)+ / dy 



1 . ^p(y) 



/o 2/ 
From Eqs. (fTl")) and (fT5j) we obtain 
^2 (a;) = h 3 (x) = ... = hg (x) , 
hi (x) + (q-l) h p (x) = 1, 

hi jx) 
Zip (a;) 

JS{(l/x) /„" d I /[/ii(»)-hp(»)]+/ !S 1 <i»[/il(»)-fcp(v)]/v}+/3H(jB) 



(15) 



(16) 



for h p (x). 

The order parameter for the Potts model is defined 
as mi = q(S(si, 1) — l/q)T/(q — 1), where (-)y is the 
thermodynamic average. So, 



&i (a:) = -[l + (<7-l)m(x)], 
9 



ft-p (x) = - [1 — m (x)] 



(17) 



These definitions satisfy the first and the second equa- 
tions of Eqs. (Tl75|) . Substituting Eqs. (fT?) into the third 
equation of Eqs. (| 16|> results in the following relation for 
the order parameter: 



m (x) 



exp 



exp 



i( - [ X dym(y) + [ dy^) illU) 

x Jo Jx y 



3[ - [dym(y) + [dy^)+ iUu) 

x Jo Jx y 



-H (18) 



V. EXACT FREE ENERGY 

For simplicity, let the external field be homogeneous. 
The field contribution to the Potts Hamiltonian © 



may be written as ~Ep=i ELo-^p^^'P) with H p — 
H5(p, 1). Let Ht be the Hamiltonian of the system at 
time t. Similarly to the Ising model, see Refs. [271. 128ll29| . 
one may obtain a recursive relation for TL t . This relation 
leads to the following recursive relation for the free en- 
ergy: 



-)3F t [H p 



E 



-0F t -ilH P +S{ St ,p)/t]+f3 £ H p 6(a t ,p) 



(19) 



»t={l,2,...?} 

So, at large t, one gets 

e -f3F t (H p )/t = e PT,UxlHv-{l/t){dF t /dH p ))5{s t ,p)_ 
s t ={l,2,...q} 

(20) 

Using hp\H p ) = dxh p (x) = —(dF t /dH p )/t gives the 
asymptotically exact free energy: 



F = -- In exp [P (h® + HS (p, 1) 



P =i 



(21) 



On the other hand, this expression may be easily ob- 
tained in the framework of the mean-field theory by sub- 
stitution Eqs. (Tnj into formula (fT5|) . which demonstrates 
the exactness of the mean-field theory. 



VI. SOLUTION OF EQ. (Tl8j) 

Self-consistent equation (fTg)) is a direct generalization 
that for the Ising model, so its solution is straightforward. 
One may see that close to x = 0, 



m(x~0) = 1 - qA(P,H)x 2f3 , 



(22) 



where A(P,H) is independent of x = i/t. If, however, 
H = 0, then m(x, T>T C ) = 0. At x = 1 (i = t), we have 



m (1) 



S /3[M+H(1)] _ 1 



e P[M+H(l)] + (g _ 1) ' 

It is convenient to define 



(23) 



n(z) = p(- [ X dym(y)+ [ dy^) + pH (x) , 

\ x Jo Jx y J 

(24) 

where z = — In x. Therefore the order parameter in terms 
of n(z) is 



m (x) 



(25) 



e «(2) + (q_ 1)' 
If H = 0, we obtain a second order differential equation: 



dn (z) d 2 n (z) 
dz dz 2 



= P 



e"( z ) + (q-l)' 



(26) 
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The boundary conditions for this equation can be derived 
from relations and ([22]) : 



( dn ( 2 



n(z = 0) = pM, 



V dz J z=0 
n(z — ► 00) = /3z + const 



(27) 



With variables n, w(n)=(3 1 dn/dz, Eq. (|2"6")l is reduced 
to a first order differential equation: 



dw fl-i f 

10-r- = p U> 



e" - 1 



(28) 



*~ V e n + (g - 1)J ' 
with the boundary conditions 

w(n = (3M) = M, 

w (n 00) = 1. (29) 

At small n, Eq. (|28p takes the form 

dui -i / 71 \ , 

"s^ 1 !*- i)- (3o) 

The physical solution of Eq. (|3"0)) satisfies 
- In [n 2 — ngw (n) + /3gu> 2 (n)] + 



1 / 2pw(n)-n \ „ , 01 , 

arctan | v y . — ) = C, (31) 



nV(W?)-l 



which gives the critical temperature /3 C = g/4. Only for 
(3 > q/4 there exists a physical solution. Near T c the 
solution of Eq. ([30]) is: 



(n,P c = l) = j[l-f(n)], 



where / (n) satisfies the relation 

1 



ln[n/(n)] + 



/(«) 



= ln< 



(32) 



(33) 



Here the constant c = 1.554... ensures that w c (n) fits 
the corresponding solution of Eq. (|28p which approaches 
1 as n — * 00. To obtain the critical singularity of the 
order parameter M(/3), one should, first, substitute the 
boundary condition w(n — (3M) — M into relation (f3Tj) . 
This gives 



C = In (/3M ) + 



1 



: arctan 



i.e., near (3 = g/4, 
C = In 



y/(40/q) - 1 ' -yy/W/q)-!)' 

(34) 



tt/2 



vW7«) - 1 



- 1. 



(35) 
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FIG. 3: The local linear response \(x,z) — dm(x) / dH (z) 
to a field applied at z = 0.5 for the 2- and 3-state Potts 
models. The values of q = 2, 3 are indicated on curves. The 
homogeneous component of the applied field is 0. The inverse 
temperature p = 1/T = 0.24 > 1/T c (q = 2,3). 



Second, near = q/A, in the region /3M <C n <C 1, 
Eq. (|3"Tj) takes the form: 



In 




+ 


' qw(n)~ 








2n 



tt/2 



= c. 

(36) 

Sewing together, from relations (|32[) . ([3^]) and (|5rj|) we 
have 



In (c) 



tt/2 



(37) 



and so, using Eqs. (|35|) and (|37)) . we readily arrive at the 
result Igj) for M(0), 

Substituting the order parameter from relation (JT7J) 
into Eq. (|21[) we obtain the free energy in terms of the 
order parameter M: 



t t 
F= In 

q P 



e «M)M/ ! + ( 9 _ 1)e - 



■0M/q 



(38) 



This leads to the result for the specific heat C(T) = 
— (T/t)d 2 F/d 2 T. Similarly, one can derive the suscepti- 
bility in zero field, Eq. ([5]). 

One may easily find a linear response of the order 
parameter at point 1 to a field applied at point z: 
x(x, z) — dm(x)/dH(z). Differentiating of (fTgj) by H (z) 
at T > T c we arrive at the following relation for the linear 
local susceptibility: 



X{x,z) = - 



- f dyx(y,z) + f —x(y,z 
x Jo Jx y 



) + S(x-z) 



(39) 
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Iterating this equation gives 



X (x, z) = x (x, z) + - S (x - z) 
(see Fig. [3] for q = 2, 3, and z = 0.5), where 



X(x < z,z) 



02 x (-l+y/l-4/3/q)/2 



i is, Aft zC-i+V 1 - 4 / 3 /?)/ 2 

,(-1-^/1-4/3/? )/2 _ *P z . 



x(x > z,z) 



p2 z (-l+y/l-4/3/q)/2 



„(-l-Vl-4/3/ 9 )/2 _ 
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4/3 jl-i+Vi^l/a 



(l + ^1- 4/3/?)' 



(40) 



(41) 



(compare with Refs. [23,|28]). So the total linear response 
to a local field applied at z, x( z ) = Jo ^ a; x( a ') z )j 1S: 

X (z) = f /q g -ii-Ji=m/2. (42) 
l + y/1- 40/q 

These formulas directly lead to results ([7]) and ((SJ in 
Sec. IIIII for the distributions of the "correlation volume" 
and of Xij- 



the traditional mean-field theory. The fact that a model 
with a first order phase transition in the usual mean-field 
regime may show the infinite order phase phase transi- 
tion with BKT critical singularity, if it is placed on the 
specific asymmetric annealed network, is the main result 
of this paper. 

Above T c (q), the correlation volume distribution P(n) 
is power-law, i.e., critical. Thus T c (q) separates the "crit- 
ical phase" and the phase [in this case T > T c (q)] where 
this distribution rapidly decreases. This is a generic fea- 
ture of the BKT transitions. Interestingly, we have found 
that the distribution P(v) of the pairwise correlations — 
the distribution of the Xij value — asymptotically coin- 
cides with P(n) above T C) which indicates that large 
Xi = J2j Xij responses ( "correlation volumes" ) are de- 
termined only by large Xij contributions. A power-law 
P{v) distribution was earlier found in the entire normal 
phase of a different cooperative model on a quite different 
network substrate, see Ref. 32]. 

Our results, generalizing the earlier findings for the 
Ising model, were obtained for a very specific annealed 
network. We suggest that resembling critical features 
may be found in the g-state Potts model on the cor- 
responding growing networks with "quenched" disorder 
and an exponentially decreasing degree distribution. Our 
findings show that the BKT critical singularity should 
occur in a wider range of cooperative systems that one 
might expect. 



VII. DISCUSSION AND SUMMARY 
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